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ABSTRACT 

We report on the X-ray spectral variability of the Seyfert 1 galaxy NGC 405 1 observed 
with the Rossi X-ray Timing Explorer (RXTE) during a 1000 day period between May 1996 
and March 1999. The spectra were obtained as part of monitoring observations and from two 
long observations using the RXTE Proportional Counter Array (PCA). During the monitor- 
ing period the 2-10 keV flux of NGC 4051 varied between lO^^^ and 7 • 10""erg/(cm2 s). 
We re-analysed RXTE PCA observations from a distinct low state in May 1998 using the 
latest background and detector response models. The RXTE and BeppoSAX observations of 
NGC 405 1 during the low state show a very hard spectrum with a strong unresolved fluores- 
cence line. This emission, probably due to reflection from a molecular torus, is likely to be 
constant over long time-scales and is therefore assumed as an underlying component at all 
flux states. By subtracting the torus component we are able to determine the spectral variabil- 
ity of the primary continuum. In the variable component we observe a strong anti-correlation 
of X-ray flux and spectral hardness in the PCA energy band. We show that the changes in 
hardness are caused by slope variability of the primary power law spectrum rather than by 
changing reflection or variable photoelectric absorption. The primary spectral index varies 
between T = 1.6 for the faintest states and F — 2.3 during the brightest states, at which level 
the spectral index approaches an asympotic value. We find that the response of the flux of the 
6.4 keV iron fluorescence line to changes in the continuum flux depends on the timescale of 
the observation. The profile of the line is very broad and indicates an origin in the innermost 
regions of the accretion disk. 
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1 INTRODUCTION 

It is thought that the central engines of active galactic nuclei (AGN) 
are powered by the infall of matter onto a supermassive black 
hole. Theoretical arguments suggest the formation of an accretion 
disk which extends from the innermost stable orbit to an outer ra- 
dius on the scale of light days. The popular unified scheme for 
Seyfert galaxies (eg Antonucci 1993) requires a dusty torus, prob- 
ably coplanar with the accretion disk, which is opaque to radia- 
tion from the near infrared to soft X-rays. According to the uni- 
fied scheme Seyfert type 1 and Seyfert type 2 AGN are intrinsi- 
cally identical object classes with the only difference being that in 
Seyfert type 2 galaxies, due to their different orientation relative 
to the observer, the central source and the broad line region (BLR) 
are obscured by a dusty torus. The first probable direct detection 
of a torus was shown in our earlier SAX (Guainazzi et al. 1998) 
and RXTE (Uttley et al. 1999) observations of NGC 4051. How- 
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ever, in general, detection of the torus and the determination of its 
geometry has not been possible. 

X-ray observations have shed some light on the structure of 
the innermost regions of AGN. The soft X-ray spectrum often has 
a relatively steep slope, which has been explained by thermal emis- 
sion from an accretion disk. However, there are alternative theoret- 
ical interpretations and due to the limited spectral resolution in this 
energy range the soft X-ray spectra are poorly understood. 

At higher energies, a hard power law spectrum dominates the 
emission. This component is probably due to Comptonization of 
the thermal UV photons in a hot corona surrounding the disk. 

The averaged Ginga medium energy X-ray spectrum of sev- 
eral Seyfert galaxies showed the presence of an emission line at 6.4 
keV, an absorption feature at 7-8 keV, and a further flattening of the 
spectrum beyond ~ 10 keV. These features are interpreted as iron 
Ka emission, iron K edge absorption and Compton reflection and 
therefore are evidence for reprocessing of X-rays by relatively cold 
matter. The discovery of broadening due to gravitational redshift 
and Doppler shifts of the iron fluorescence lines in Seyfert X- 
ray spectra suggests that at least part of the reprocessing takes place 
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in the inner parts of the accretion disk ( Tanaka et al. 1995, Nandra 
et al. 1997, Reynolds & Begelman 1997). However, the detection 
of narrow components to the Fe Ka Une with ASCA (eg Weaver et 
al. 1997) and recently Chandra (Yaqoob et al. 2001, Kaspi et al. 
2002) and XMM-Newton (Reeves et al. 2001) suggests reprocess- 
ing of X-rays in matter at a larger distance from the black hole, eg 
in the BLR or in the molecular torus. 

The analysis of X-ray variabiUty is a powerful tool for the in- 
vestigation of the inner regions of the AGN central engine, since the 
time-scales of the variability can give an indication of the geomet- 
rical sizes of the regions involved. Yaqoob et al. (1996) reported 
a rapid (< 3 • 10* s) response of the Fe fluorescence line flux to 
the continuum flux in NGC 7314. From the resulting maximum 
distance of the reflecting material from the primary X-ray source 
they derive an upper limit on the mass of the black hole. However, 
tdhe analysis of fluorescence line variability in a number of other 
Seyfert galaxies have led to confusing and partly contradictory re- 
sults. While in some sources a positive correlation of line and con- 
tinuum flux has been observed, averaged over long timescales, (eg 
NGC 5506, Lamer, Uttley & M'^Hardy 2000), other studies found 
no evidence for a close relationship between continuum and fluo- 
rescence line (Weaver, Gelbord & Yaqoob 2001, Vaughan & Edel- 
son 2001). 

Flux-dependent variability of the continuum spectral shape it- 
self can be used to constrain models of the Comptonising accretion 
disk corona (eg Haardt, Maraschi & Ghisellini 1997). In general, 
the continuum slope appears to be correlated with flux (eg Leighly 
et al. 1996; M'^Hardy, Papadakis & Uttley 1998 ; Lamer, Uttley 
& M'^Hardy 2000), as predicted by most simple Comptonisation 
models. However, it is important to determine whether the degree 
of slope variability, and the form of the correlation with flux (eg 
does the continuum slope saturate at some maximum value?) are in 
agreement with existing models. 

NGC 4051 is a nearby (z=0.0023) low luminosity Seyfert 1 
galaxy, which is among the most variable AGN in the X-ray band 
(Green, M'^Hardy & Lehto 1993). We have been monitoring NGC 
4051 with RXTE since 1996. Since our monitoring commenced, 
the 2-10 keV flux from the object has varied by a factor of ~100. 
In 1998 NGC 4051 entered a state of extremely low X-ray flux that 
lasted for ~ 150 days (Guainazzi et al. 1998, Uttley et al. 1999). 
Here we report on the spectral changes that accompanied the dra- 
matic flux variability. In section 2 we describe the observations and 
the reduction of the RXTE data. Results from observations in May 
1998, during the extreme low state, are re-analysed in a consistent 
manner with the other observations and are discussed in section 3. 
In section 4 we present the analysis of the continuum and fluores- 
cence line variations observed during the long term monitoring and 
during an RXTE long look in December 1996. 



2 OBSERVATIONS AND DATA REDUCTION 

We have been monitoring NGC 4051 with RXTE since May 1996. 
The observations were separated by a range of time intervals. For 
the first 6 months the observations took place weekly with periods 
of twice-daily observations in May 1996 and daily observations in 
October 1996. Since November 1996 the source has been observed 
at fortnightly intervals. RXTE observed NGC 4051 wifli tiie Pro- 
portional Counter Array (PCA) and the High Energy X-ray Timing 
Experiment (HEXTE). The PCA (Zhang et al. 1993) consists of 
5 Xenon-filled Proportional Counter Units (PCUs), sensitive to X- 
ray energies from 2-60 keV. The maximum effective area of the 



PCA is 6500 cm'^. As the signal to noise ratios of the higher energy 
HEXTE spectra on NGC 405 1 are low, we only present the results 
from the PCA. For technical reasons the high voltage settings of 
the PCA units have been adjusted several times during their hfe- 
timc. In this paper we only present data taken during the PCA gain 
epoch 3, which lasted from April 1996 to March 1999. The total ex- 
posure time of the 130 individual monitoring observations during 
gain epoch 3 is 85 ksec. Apart from the monitoring observations 
RXTE performed two long-look observations of NGC 4051 in De- 
cember 1996 with a total on-source time of 36.4 ksec within 3 days 
and in May 1998 for 84.4 ksec (on-source) within 2 days. 

We have used ftools v4.2 for the reduction of the PCA and 
HEXTE data. PCA "good times" have been selected from the Stan- 
dard 2 mode data sets using the following criteria: target elevation 
> 10°, pointing offset < 0.01°, time since SAA passage > 30 min, 
standard threshold for electron contamination. Occasionally one or 
more of the 5 Proportional Counter Units (PCUs) of the PCA were 
switched off during the observations. In order to maximise the sig- 
nal to noise ratio of the spectra, we used the data from all PCUs 
that were switched on during each individual pointing and used 
only PCU layer 1. We calculated the background in the PCA with 
the tool PCABACKEST v2.1 using the L7 model for faint sources, 
which is suitable for determining the PCA background for energies 
< 24 keV. From a series of 100 blank field observations from the 
public archive (proposal ID P30801) we determined the uncertain- 
ties in the estimated background rates as a function of energy. At 
energies below 24 keV the Icr uncertainty of the background is of 
the order 0. Diets s^^kcV"^ (see Lamer, Uttley & M'^Hardy 2000 
for details). The PCA response matrices were calculated invididu- 
ally for each observation using pcarsp v2.37, taking into accoimt 
temporal variation of the detector gain and the changing numbers 
of detectors used. 



3 THE LOW STATE SPECTRUM 

Our monitoring with the RXTE PCA revealed that during a 150 day 
period from January 1998 to May 1998 NGC 4051 was in a state 
of imusually low X-ray flux (Uttley et al. 1999, see Fig. 1). Near 
the end of this period NGC 4051 was observed simultaneously by 
RXTE (Uttley et al. 1999) and BeppoSAX Guainazzi et al. (1998). 
We derived a 2-10 keV flux of 2.0 • 10"^^ergs"^cm-^ from the 
RXTE PCA spectrum and 1.4 • 10"^^ergs"^cm"^ from the Bep- 
poSAX MECS data. Both observations showed a very hard X-ray 
spectrum and unresolved 6.4 keV iron line emission with a very 
high equivalent width (~ 1000 eV). This spectrum was interpreted 
as a pure reflection component with the central source being vir- 
tually switched off. Since the reflected emission was still present 
after the source had been in a low state for 150 days, the distance 
of the reflector from the central source must be > 10^^ cm. It is 
likely that the Compton reflection and iron line fluorescence from 
the distant reflector is a constant contribution to the X-ray spectrum 
of NGC 405 1 . For the investigation of the variabiUty of the primary 
X-ray source in NGC 405 1 it is necessary to take this underlying 
component into account. 

We have therefore re-analysed the PCA data from the low state 
using the latest response matrices, appropriate to gain epoch 3, and 
background models. We fitted two models to the PCA and Bep- 
poSAX MECS spectra: a simple power law model with a Gaus- 
sian emission line and a pure reflection model, also with a Gaus- 
sian emission line. The model for a pure reflection component was 
calculated by using the XSPEC PEXRAV model (Magdziarz & 
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Figure 1. 2-10 keV Ughtcurve of NGC 4051 during the long term monitoring campaign (filled circles with error bars). The diamonds show the lightcurve from 
the December 1996 observation, the open circle indicates the flux level during the May 1998 observation at the end of the 150 day long low state. 



Table 1. Fit parameters for the May 1998 low state spectrum 



model 


r 




f6.4 


EW6.4 


EW7.057 


EW7.477 


(dof) 








IQ-^ph 


[keV] 


[keV] 


[keV] 










RXTE PCA 










PWL+ZGAUSS 


0.56 it 


11 


2.80 ± 0.70 


1.40 ±0.35 






38.8 (42) 


PWL+3 ZGAUSS 


0.52 ± 


11 


2.57 ±0.53 


1.40 ± 0.34 


0.16 


0.10 


36.8 (42) 


PEXRAV+ZGAUSS 


2.02 ±0 


10 


2.16 ±0.75 


0.79 ± 0.35 






39.8 (42) 


PEXRAV+3 ZGAUSS 


2.01 ±0 


10 


2.35 ± 0.66 


0.88 ± 0.27 


0.12 


0.11 


38.0 (42) 








BeppoSAX MECS 










PWL+ZGAUSS 


0.37 ±0 


17 


1.88 ±0.60 


1.16 ±0.37 






28.4 (29) 


PWL+ 3 ZGAUSS 


0.52 ± 


18 


2.04 ±0.64 


1.26 ±0.39 


0.14 


0.09 


25.5 (29) 


PEXRAV+ZGAUSS 


1.75 ±0 


18 


1.85 ± 0.64 


0.87 ±0.30 






33.9 (29) 


PEXRAV+3 ZGAUSS 


1.83 ±0 


20 


1.12 ±0.63 


0.50 ±0.28 


0.07 


0.06 


29.6 (29) 



Zdziarski 1995) with a fixed reflected fraction R = 100. With this 
setting the contribution of the intrinsic power law emission is neg- 
ligible. See Table 1 for the results of the spectral fits. Throughout 
the paper 1 a confidence limits are given as error values. The er- 
ror of the fit parameters in Table 1 were calculated for 4 interesting 
parameters. Due to the relatively low signal-to-noise ratios of the 
data a clear cut decision between these two models is not possible. 
For both the PCA and MECS spectra the power law model yields 
a slightly better value than the pure reflection fit. However, the 
power law fit requires a photon index of F ~ 0.5, a value that is un- 
likely for standard Comptonization models but which might be ap- 
plicable to an advective flow. If the low state X-ray emission solely 
arises from reflection, the PEXRAV model fit to the PCA spectrum 
results in a photon index F ~ 2.0 for the primary emission. Assum- 
ing that the reflecting matter subtends a solid angle of Itt steradians 
or less, the required incident radiation corresponds to a 2-10 keV 
flux level of > 2 • 10~'^'^ergs~^cm~^. This is in good agreement 
with the flux range of ~ 1 ■ 10^^^ to 6 ■ lO^^^crgs^^crri^^ 
found during the monitoring of NGC 4051 before the 150 day low 
state period. 



The strong iron Ka fluorescence line supports the interpreta- 
tion that the low state spectrum is a pure reflection component. The 
very low level of the continuum in these observations of NGC 405 1 
allow a more detailed analysis of any fluorescence or emission 
lines. When fitting a power law or PEXRAV model with a single 
Gaussian line to the PCA spectrum, the resulting source frame line 
energy of ~ 6.6 ± 0.1 keV is slightly higher than the nominal value 
of 6.4 keV for fluorescence from cold iron. The best fit model to the 
BeppoSAX MECS spectrum includes a narrow line at 6.47 ± 0.06 
keV, consistent with 6.4 keV. The fit residuals, however, show some 
excess flux at ~ 7 keV. When fixing the energy of the fluorescence 
line to 6.4 keV, excess flux at ~ 7 keV is evident in all spectra. We 
therefore included narrow (0.05 keV) iron Ka and nickel fluo- 
rescence lines at 7.058 keV and 7.477 keV into the spectral model. 
The fluxes of these additional lines were set to Fe K/3/Fe Ka = 
0. 1 1 as theoretically expected for neutral iron (Kikoin 1976) and 
Ni Kq/Fc Kq, = 0.07 as expected for solar abundances (Nandra et 
al. 1999). Taking these lines into account slightly improves the 
of both the PCA and the MECS fits without increasing the number 
of free parameters (see Table 1). 
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Figure 2. Huorescence line profile in the lowest flux state during the De- 
cember 1996 observation. The flux of the narrow component is consistent 
with the Une flux during the May 1998 low state. 



Throughout the remainder of this paper, we will assume that 
a narrow iron line component, whose flux is consistent with the 
line flux measured in May 1998, is present in the spectrum of 
NGC 405 1 over the 3 years of our RXTE observations. Fig. 2 shows 
the Fe Kq line profile with a possible strong narrow component 
during the lowest flux states of the December 1996 observation 
(see first row of Table 3). This is consistent with our view that a 
component, which is constant on the time-scale of years is present 
in the X-ray spectrum of NGC 4051. In the analysis of the spectral 
variabihty in the RXTE PCA data, we therefore add the PEXRAV 
-I- 3 ZGAUSS model with the parameters derived from the PCA 
low state data (see Table 1) as a constant component to all model 
fits of the monitoring and December 1996 data. 

We note that the presence of narrow line components in 
some other objects has been confirmed by Chandra and XMM ob- 
servations. Yaqoob et al. (2001) detected a narrow line at 6.40 
keV in a Chandra High-Energy Transmission Grating spectrum of 
NGC 5548. Since the line is slightly resolved, they assume an ori- 
gin in the broad line region. Ogle et al. (2000) suggested that part 
of the narrow iron Une observed by Chandra in NGC 4151 may 
originate in the narrow-line region. Reeves et al. (2001) found an 
unresolved component at 6.4 keV in the XMM EPIC spectrum of 
the low-luminosity quasar MKN 205. Kaspi et al. (2002) were able 
to resolve the narrow 6.4 keV Une in NGC 3783 with the Chandra 
HETGS. The velocity dispersion of FWHM=1720 kms"^ suggests 
that the Une is emitted in the outer BLR or the inner part of the 
molecular torus. 

Although we can say little about the line width from our 
RXTE observations, the fact that the line remains detectable in the 
very low state May 1998 observations, 150 days after the contin- 
uimi went into the low state, means that it could not come from the 
broad line region, but part of it might come from the narrow line 
region. However the very hard continuimi is not easily explained 
by scattering from the narrow Une region. 

If the strong torus contribution to the iron fluorescence line 
in NGC 4051 is common in other Seyfert galaxies, some of the 
previous results on disk line spectroscopy will have to be revised, 
once high resolution spectra become available. 



Table 2. Flux bins for monitoring observations 





flux range 


mean flux 


exposure 




[erg s" 1 


cm ^] 




[sec] 




< 1.0 • 10-" 


5.47 


10-12 


21824 


1.0 


-2.0-10"" 


1.63 


10-11 


22352 


2.0 


-3.0 -10-" 


2.76 


10-11 


15744 


3.0 


-3.5 -10"" 


3.36 


10-11 


8288 


3.5 


-4.5 • 10"" 


4.06 


10-" 


10880 


4.5 


- 5.5 ■ 10-" 


5.21 


10-11 


4176 




> 5.5 ■ 10-" 


6.66 


10-11 


2016 



Table 3. Flux bins for December 1996 observation 

flux range mean flux exposure 



[erg s 


-1 


cm 2] 






[sec] 


< 1.0 ■ 10- 


11 


6.78 


10- 


12 


20800 


1.0 - 2.0 ■ 10- 


11 


1.53 


10- 


11 


7792 


2.0 - 3.0 • 10- 


11 


2.65 


10- 


11 


6800 


> 3.0-10- 


11 


3.62 


10" 


11 


1008 



4 SPECTRAL VARIABILITY 

The large variations in the X-ray flux observed within the 3 years 
of monitoring of NGC 4051 make the RXTE observations ideally 
suited for an investigation of associated spectral variability. Since 
the signal to noise ratio in any one of the individual 1 ksec pointings 
is not sufficient for individual spectral analysis, we have grouped 
and added the spectra according to source flux. Since the data were 
gathered by different combinations of PCUs and the detector gains 
changed slightly over time, we used fitted model fluxes for the 
grouping of the data. In order to obtain the flux level for each spec- 
trum we fitted a single power law model to each spectrum in the 
2-10 keV range using the appropriate detector response for each 
spectrum. The monitoring observations were grouped into 7 flux 
bins, while the spectra from the December 1996 long look obser- 
vation were grouped into 4 bins (see Tables 2 and 3 for the bound- 
aries of the bins). The appropriate detector response matrices for 
each spectnmi have been created by averaging the response matri- 
ces of the individual input spectra (using their integration times as 
weighting factors). 

As described in section 3 we regard the emission observed 
during the low state in May 1998 as a constant underlying spectral 
component from a distant reflector. For all spectral fits in this sec- 
tion we have therefore included the best fit of the PEXRAV -i- 3 
ZGAUSS model to the May 1998 PCA data as a fixed additional 
component. 

4.1 Continuum and reflection 

Simple power law fits to the flux binned spectra from the moni- 
toring observations after removal of the torus component still re- 
veal a strong correlation of the spectral slope with the flux level of 
NGC 4051: the softest spectra are found when the flux is highest 
(see fig 7). In the context of the standard models for the X-ray emis- 
sion of Seyfert galaxies three principal mechanisms can contribute 
to this correlation: 

1. Intrinsic spectral variability of the primary continuum. During 
states of high luminosity the higher flux of UV/soft X-ray photons 
may cool the accretion disk corona and therefore lead to a softer 
Compton spectnmi at higher X-ray energies. 
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2. The fraction of reflected radiation may be higher during peri- 
ods of low luminosity (eg due to an additional, non- variable reflec- 
tion component). 

3. Gas near the X-ray source that is usually strongly ionised by 
radiation from the central source might recombine to less ionised 
states when the luminosity of the AGN drops. The resulting in- 
crease in absorption of soft X-ray then leads to a hardening of the 
spectrum. Outside of the very low luminosity states NGC 4051 is 
observed to have a signiflcant "warm absorber" with a column den- 
sity of 5 - 8 • lO^^crrr^ ( M'^Hardy et al. 1995, Komossa & Fink 
1999) which could lead to substantial absorption of soft X-rays if 
the warm absorbing gas recombined. 

In order to distinguish between the flrst two possibilities we 
fitted PEXRAV models with a variable Gaussian line in addition 
to the fixed torus reflection model to the flux selected spectra of 
both the monitoring and December 1996 observations. We reduced 
the number of free parameters by fixing some of the less critical 
parameters of the PEXRAV model to reasonable values. The high 
energy cutoff of the primary power law was set to 300 keV, effec- 
tively eliminating the cutoff. The iron and metal abundances were 
set to the solar values and the inclination angle was fixed at 30° 
as suggested by the results of disk line model fits to ASCA spectra 
(Nandra et al. 1997). Confidence contours of the primary photon 
index F and the reflected fraction R are shown in Figs 3 and 4. It 
is evident that a continuum of variable spectral index is needed to 
satisfy the data. The best fit photon indices range from F = 1.6 for 
the lowest flux states to F = 2.3 for the highest flux states of the 
monitoring observations. Note that the inclusion of the presumed 
constant torus component does not change the results qualitatively. 
If we omit this fixed component, the photon index of the lowest 
flux spectrum is F = 1.2; the fit to the highest flux spectrum is 
not affected by this change. No variabiUty of the reflected fraction 
is obvious and at the 90% confidence levels all spectra are consi- 
sistent with zero reflection. Apart from the fact that the constant 
reflected (torus) component will become less prominent at higher 
flux levels, there is no obvious contribution of reflection to the spec- 
tral variability. 

We also investigated the warm absorber as the third possible 
source of hardness variability. Again including the low state spec- 
trum as a constant component, we modelled the flux-selected X-ray 
spectra with a power law and Gaussian emission line absorbed by a 
warm absorber (XSPEC ABSORI model) with a column density 
of A^H = 5 • lO^^cm"^, as suggested by Komossa & Fink (1999). 
The models were calculated on a grid in spectral index F and ioni- 
sation parameter ^ and fitted to each of the 7 flux selected spectra of 
the monitoring observations. As a result, the ionisation parameter 
appears to be poorly constrained. The spectral index is clearly cor- 
related with flux and there is no F that would be consistent with all 
of the 7 spectra (see Table 4). We therefore conclude that changes 
in the ionisation state of the warm absorber in NGC 405 1 do not 
play a signiflcant role in the spectral variability of the source. 

To summarize, variability of primary power law continuum 
slope is the only viable explanation of the strong hardness variabil- 
ity in the 2-24 keV X-ray spectrum of NGC 4051. 

4.2 Fluorescence line variability 

After subtraction of the low state reflection spectrum with the nar- 
row iron fluorescence line, broad fluorescence line emission is still 
detectable at all flux levels. Figs. 5 and 6 show ratio plots obtained 
by fitting power law models to the flux selected spectra excluding 



Table 4. Results of warm absorber fits to flux selected monitoring spectra 



mean flux 




photon index 




[ergs 


^ cm 2] 


la range 


1(7 limit 


5.47 


10- 


12 


1.50-1.91 


>0 


1.63 


10- 


11 


1.85-2.14 


> 


2.76 


10- 


11 


2.10-2.30 


> 60 


3.36 


10- 


11 


2.13-2.28 


> 120 


4.06 


10- 


11 


2.08-2.28 


> 370 


5.21 


10- 


11 


2.28-2.43 


> 170 


6.66 


10- 


11 


2.29-2.45 


> 240 




1.0 1.5 2.0 2.5 3.0 

r 



Figure 3. 68% and 90% confidence contours of the parameters photon in- 
dex r and reflected fraction R for the PEXRAV model fits to the flux se- 
lected PCA spectra. The letters 'a' to 'g' denote the spectra in the order of 
increasing flux levels. See section 4 for a fuU description of the model fits. 

the 4-8 keV spectral range. The line is broad and strongly redshifted 
as expected for an origin in the inner parts of an accretion disk. 

We fitted each of the flux selected spectra with a model con- 
sisting of the best fit low state spectrum as discussed in section 3, a 
power law, and the XSPEC DISKLINE model. The DISKLINE 
model calculates the profile of a line originating from a Keplerian 
disk surrounding a Schwarzschild black hole. The rest energy of 
the line was fixed at 6.4 keV. The geometric origin of the fluores- 
cence line is defined by three parameters: the inner radius -Ri, the 
outer radius Ro and the power law index q which describes the ra- 



3 




1.0 1.5 2.0 2.5 3.0 



Figure 4. Same as Fig. 3 for the flux selected spectra of the December 1996 
observation. The letters 'a' to 'd' denote the 4 flux levels in ascending order. 
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Figure 5. Iron line profiles in the flux selected spectra of the monitoring 
observations. The letters a-g denote the 7 flux intervals in ascending order. 
The low state spectrum including the narrow 6.4 keV fluorescence line has 
been subtracted. Note the strong broadening of the line at aU flux levels. 



dial variation of the emissivity. We set Ri to 6 rg, the radius of the 
innermost stable orbit for a Schwarzschild black hole. The outer 
radius Ro was set to 400 rg. The emissivity index q was left free to 
vary. 

Table 5 summarizes the results of the disk line modelling. For 
each spectnmi we give both the equivalent widths of the narrow flu- 
orescence lines from the fixed narrow line component and the best 
fit values for the broad disk line. The values for the narrow lines are 
derived from the total fluxes of the lines at 6.4 keV, 7.06 keV, and 
7.477 keV (see section 3). Note that the narrow line fluorescence 
from the torus can be a major contribution to the total line flux, in 
particular during lower flux states. Hence the profile of the disk line 
can only be determined if the contribution from the torus can be re- 
liably measured. The best fit values of the disk inclination confirm 
the values of i ~ 30° resulting from ASCA observations (Nandra 
et al. 1997). The values of q in most cases are poorly constrained. 
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Figure 6. Iron line profiles for the December 1996 observations in the 4 
different flux ranges. As in Figure 5 the low state spectrum including the 
narrow line has been subtracted. 



but generally the emissivity drops steeper than with q = —3. This 
indicates a concentration of the line emission in the innermost re- 
gions of the disk, as suggested by the broad profile of the iron line. 
Since the PEXRAV model fits described in section 4 required no 
disk reflection (see Figs 3-1-4) we do not include disk reflection in 
the disk line model fits. However, we repeated the model fits using 
a PEXRAV model with reflected fraction R = 1 instead of the 
power law model. Including reflection to the fits did not change the 
line fluxes or disk line parameters significantly. 

The best fit photon indices, line fluxes, and line equiva- 
lent widths are plotted in Fig. 7. The best fit values of the 
fluorescence line flux from the long term monitoring observa- 
tions show a clear correlation with the 2-10 keV flux. However, 
the correlation is not linear, in the continuum flux range range 
(2.()...5.5)10^^^crg cm~^s~^ the line flux is nearly constant. The 
variation of the equivalent width is not very pronounced, the lowest 
and the three highest flux bins are consistent with EW = 200eV. 

On the shorter time scales of the December 1996 observation 
the Une equivalent width decreases with continuimi flux, although 
the line flux and continuimi flux are still correlated. 



5 DISCUSSION 

We show that the variable X-ray spectrum of NGC 4051 can be 
well described by a model that includes two principal components: 

1. A hard, constant, component including a narrow iron fluores- 
cence line as revealed during the extreme low state in May 1998. 

2. A variable component including a power law with strongly 
variable slope and a very broad emission line. We find a strong cor- 
relation between the power law slope and the source flux. On the 
long time-scales of the monitoring observations the Une flux is cor- 
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Table 5. Disk line model fits to flux selected spectra 
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Figure 7. Flux dependencies of power law spectral index (top), disk line 
flux (middle) and disk line equivalent width (bottom). The results from the 
monitoring observations are plotted as dots, the December 1996 long look 
is represented by open squares. Error bars indicate Ax^ = 1.0. 



related with the continuum flux, although not in a simple manner. 
The reflected fraction of the disk component is less than i? = 1 at 
all flux levels. 

This two-component model is consistent with the model we 
proposed for the X-ray spectral variability of the Seyfert 2 galaxy 
NGC 5506 (Lamer, Uttley & M'^Hardy 2000). 

We note that even if the constant hard component is not in- 
cluded in our spectral fits, the general form of our results is not sig- 
nificantly affected, as the assumed constant component is relatively 



weak. Indeed, the lack of inclusion of a hard continuum component 
would result in an even larger spectral index variation with flux. 



5.1 The Iron Line 

Even on the long timescales probed by our time-averaged moni- 
toring data, the flux-dependent behaviour of the broad iron line is 
complex. For these data, the line flux increases more-or-less pro- 
portionally with the continuum flux, resulting in a roughly con- 
stant equivalent width over a decade range of flux. The Decem- 
ber 1996 long-look data also show line flux increasing with con- 
tinuum flux, although the relation is not directly proportionate, so 
that the equivalent width decreases with flux. In fact, the line fluxes 
measured in December 1996 are consistently larger than the corre- 
sponding fluxes measured from the long-term monitoring data. The 
anti-correlation of Une equivalent width and continuum flux in De- 
cember 1996 timescales is also in contrast to the result of (Wang et 
al. 1999), who report a positive correlation during an A5CA obser- 
vation in 1994. 

The discrepancy between the iron line behavioiu' in the De- 
cember 1996 and long-term monitoring data might be explained if 
there is additional short-term variability in the iron line which is 
not simply related to the continuum flux. For example, Vaughan & 
Edelson (2001) show fliat in flie Seyfert 1 MCG-6-30-15, the broad 
iron line flux varies significantly but independently of short term 
continuum variations. 

One possibility is that the iron line flux tracks the long-term 
variations in the continuum flux (which are being probed to some 
extent with the long-term monitoring data), but responds only 
weakly to the short-term variations which are observed during the 
December 1996 long-look observation. A number of theoretical pa- 
pers have been written to explain why the iron line flux may not 
vary linearly with the continuum flux (eg Matt et al. 1993, Nayak- 
shin & Kazanas 2002, Ballantyne & Ross 2002), often involving 
ionised discs but these models have so far been largely untroubled 
by data. We are aquiring more long-term monitoring data, sampling 
a broader range of long-term flux variations, to determine whether 
the iron line does follow the continuum on long timescales and to 
provide some constraints for theoretical models. 
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5.2 The Reflected Component 

Our observations of NGC 405 1 do not support the correlation be- 
tween the photon index and the reflected fraction R in the PEXRAV 
model as reported from ASCA spectroscopy of a sample of Seyfert 
galaxies (Zdziarski, Lubinski & Smith 1999) . From Fig. 4 it is 
obvious that the reflected fraction remains below R = 1 even for 
the softest states of the source. There is also no evidence for this 
correlation in the RXTE spectra of NGC 5506 (Lamer, Uttley & 
M'^ Hardy 2000). We therefore suggest that the reported correlation 
does not apply to the variations of photon index and reflected frac- 
tion in a given object. 



5.3 Continuum Spectral Variability 

During our monitoring campaign the primary continuum 2-10 keV 
photon index, F, varies strongly with photon flux from 1.60 at the 
lowest flux levels to 2.35 at the highest fluxes. The same correla- 
tion is also observed on the shorter time-scales of the December 
1996 long look. Flux-F correlations have been observed before in 
NGC 4051 (Matsuoka et al. 1987) and other Seyfert galaxies (eg 
NGC 4151, Perola et al. 1986) but, as in the present paper, it has 
only recently been possible to disentangle the effects of reflection 
and variations of the primary X-ray spectrum (eg Lamer, Uttley & 
M'^Hardy 2000, Chiang et al. 2000, Lee et al. 2000). 

The slope-luminosity correlation is often explained by 
stronger cooling of the accretion disk corona during episodes 
of high thermal seed photon flux from the accretion disk itself 
(eg Pietrini & Krolik 1995, Malzac & Jourdain 2000). Haardt, 
Maraschi & Ghisellini (1997) have calculated luminosity - spec- 
tral index relations in Compton cooled accretion disk coronae. For 
a compact, pair dominated corona they predict spectral index vari- 
ations of AF ~ 0.3 for luminosity variations by more than a factor 
of 20. However the variations seen here exceed their predictions 
and imply, in their scenario, a non-pair dominated corona. In cer- 
tain regimes this model predicts a positve correlation of 2-10 keV 
flux and spectral hardness, which is not observed in NGC 4051. 
Pietrini & Krolik (1995) point out that the spectral index depends 
almost solely on the ratio of seed photon compactness Is and 
hot plasma heating rate compactness Ih with a = 1.6{la/lh)^^*. 
The observed spectral indices in NGC 4051 then correspond to 
h/lh = 0.02..0.4. 

Examination of fig 7 shows that the change of spectral index 
with flux is not linear. The rate of increase of spectral index with 
flux is very rapid at low fluxes but decreases at the highest fluxes 
where the spectral index approaches an asymptotic level. This sat- 
uration of the 'spectral index/flux' relationship has been known 
for some time; eg the saturation was clearly visible in our early 
RXTE monitoring observations of MCG-6-30-15 and was reported 
by M'^Hardy, Papadakis & Uttley (1998) where it was suggested 
that the relationship might derive from the combination of a con- 
stant spectrum hard component, and a steeper spectrum variable 
component. 

Saturation of the spectral index/flux relationship was again re- 
ported in MCG-6-30-15 by Shih et al. (2002) from a long ASCA 
observation. In MCG-6-30-15 both the long term RXTE monitor- 
ing and short term ASCA observations agree that the saturation 
level of the spectral index is ^ 2.1 (see fig 7 of M'^Hardy, Pa- 
padakis & Uttley 1998 and fig 8 of Shih et al. 2002). However the 
monitoring observations cover a wider flux and spectral range and 
show variation of the spectral index between 1 .65 and 2.05 whereas 
the continuous ASCA observation only shows an index variation 



between 1.9 and 2.1. Similariy in NGC 4051 (fig 7) we see that the 
monitoring observations cover a wider flux and spectral range than 
the December 1996 long look and, as with MCG-6-30-15, the re- 
sultant time-averaged spectral index/flux relationship is smoother. 
We note, however, that although the lowest spectral index so far 
measured in the RXTE monitoring observations is about the same 
in both NGC 4051 and MCG-6-30-15, the saturation level is ~ 2.4 
in NGC 4051 compared to ~ 2.1 in MCG-6-30-15. 

In M'^Hardy, Papadakis & Uttley (1998) we suggested that 
the torus might be the source of the possible hard, constant, compo- 
nent. However the very hard spectral component found in the May 
1998 very low state, which represents an upper limit to the torus 
contribution, was removed before producing the spectral index/flux 
relationship (fig 7). Thus, if we wish to retain the two-component 
spectral model, we require a different location for the bulk of the 
hard component. However any hard component produced closer to 
the central source by reprocessing of primary radiation is likely 
to vary, although probably not as rapidly as the primary contin- 
uum. Our time-averaged spectra wash out short term variability so 
a hard component varying slowly in line with the primary contin- 
uum could not produce the observed spectral index flux relation- 
ship. Unless a hard component is produced by some mechanism 
other than reprocessing it is therefore more likely that the change 
of spectral index with flux is driven by some change in the physical 
properties of the primary continuum source, eg by changes in the 
seed photon populations as discussed above. 

Any such physical model of the primary continuum must be 
able to produce different saturation levels in different sources and 
must be able to account for both the wide spectral range, and the 
very hard spectra measured at the lowest flux levels, in our moni- 
toring observations. 
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